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1.  Introduction 


The  analysis  of  the  radiating  infinite  flange  problem  has  been  widely  studied  and  different 
methods  for  determining  the  aperture  fields  are  presented  in  the  published  literature  (7-3).  In 
practice,  flush  mounted  aperture  antennas  are  widely  used  and  are  often  approximated  by  an 
aperture  in  an  infinite  conducting  surface  (infinite  flange).  Furthermore,  the  solutions  determined 
by  these  methods  reasonably  approximate  those  of  a  radiating  waveguide  with  no  infinite 
conducting  surface  (4).  Although  accurate,  these  methods  require  the  rigorous  calculation  of  the 
fields  at  the  aperture  of  the  waveguide. 

This  report  illustrates  a  more  straightforward  approach  that  directly  computes  the  radiated  far 
field  eliminating  the  direct  calculation  of  the  aperture  field.  This  approach  uses  a  modal 
decomposition  matrix  (MDM)  method  based  on  the  sampling  of  the  transverse  free  space  wave 
number  to  determine  the  spectral  coefficients  of  the  widely  used  stationary  phase  approach.  The 
results  determined  by  our  method  are  compared  to  the  far- field  radiation  patterns  achieved  using 
CST  Studio  Suite  software  for  the  same  problem. 


2.  Modal  Decomposition  of  a  Rectangular  Waveguide 


This  section  describes  the  theory  of  matching  the  transverse  electric  and  magnetic  fields  that 
exist  inside  a  uniform  rectangular  waveguide  to  those  of  the  radiated  far  fields  in  free  space.  The 
modal  decomposition  described  here  is  based  on  the  theory  published  by  Felsen  and  Marcuvitz 
(5,  9).  The  transverse  field  inside  the  guide  is  assumed  to  be  generated  from  the  dominant 
propagating  mode  inside  the  waveguide.  The  dominant  mode  is  determined  by  the  lowest  cutoff 
frequency  in  the  waveguide  and  for  a  rectangular  waveguide  the  transverse  electric  (TEio)  mode 
is  the  dominant  mode. 

By  matching  tangential  boundary  conditions  at  the  aperture,  a  system  of  equations  is  derived  that 
will  directly  yield  the  Fourier  transform  of  the  transverse  aperture  field.  Based  on  this,  we  can 
use  the  stationary  phase  approximation  to  calculate  the  far- field  radiation  patterns. 

Figure  1  shows  two  orientations  of  the  same  rectangular  waveguide.  A  half-space  boundary 
exists  at  the  aperture  (z  =  0),  as  depicted  in  figure  la,  and  the  waveguide  extends  to  -oo  in  the 
negative  z-direction.  This  is  known  as  the  semi-infinite  waveguide  approximation  (5).  Using  this 
approximation,  we  can  construct  the  form  of  the  electric  and  magnetic  fields  that  exist  inside  the 
waveguide  due  to  an  incident  TEio  mode. 
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Figure  1.  (a)  Illustration  of  the  boundary  between  the  waveguide  and 
free  space  at  z  =  0.  The  thickness  of  the  PEC  walls  is 
assumed  to  be  negligible  in  comparison  to  the  wavelength. 

(b)  Illustration  of  a  transverse  cross  section  of  the  waveguide. 
Again  the  thickness  of  the  ground  plane  around  the  waveguide 
is  assumed  to  be  negligible  in  comparison  to  wavelength, 
i.e.,  less  than  the  skin  depth. 


2.1  Incident  TE]0  Mode  Case 


The  incident  TEio  mode  shown  in  figure  la  is  assumed  to  exist  inside  the  waveguide  with  the 
following  form: 


ETinc{r)  =  e”N(x,y)V{z) 

Kmc{r)  =  h”N(x,y)l(z) 


(la) 

(lb) 


where  V(z)  and  I(z)  are  the  voltage  and  current  at  point  z  inside  the  waveguide,  and  are  defined 
by  solutions  to  the  wave  equation.  The  subscript  N  denotes  that  this  is  the  incident  mode 
impingent  on  the  aperture  (z  =  0)  and  the  superscript  ”  denotes  that  this  is  a  TE  mode.  The 
superscript  ’  is  used  to  denote  the  transverse  magnetic  (TM)  mode.  The  mode  functions  eD(x,y) 
and  hD(x,y)  are  defined  as 
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where  $u(x, y)  =  A'v  sin 
x¥u(x,y)  =  A”u  cos 
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,  and  i)  represents  the  (m,n)  pair  known  as  the 


mode  number  (5).  The  two  mode  number  indices  can  take  integer  values  zero  or  greater.  For 
instance,  the  TEio  mode  has  mode  indices  of  m  =  1  and  n  =  0.  A  and  A  ”v  are  determined  by 
normalizing  equation  2  across  the  transverse  plane  of  the  waveguide  in  figure  lb  as  (5) 
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where  d  is  defined  as  dmk  =  0  for  m  ±  k  and  Smm  =  0.  Solving  equation  3  for  the  normalization 
constants  yields 
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where  Pmn  is  the  incident  mode’s  amplitude.  Taking  equations  2  and  4  into  account,  we  can 
construct  the  total  transverse  E-  and  H-fields  inside  the  waveguide 
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where  r  =  x-x0  +  yy0  +  z-za,  the  subscript  M  denotes  all  mode  numbers  that  are  not  incident  on  the 
aperture,  and  r  is  the  reflection  coefficient  at  the  aperture.  Notice  that  both  the  incident  and  non¬ 
incident  modes  are  reflected  from  the  aperture  at  z  =  0.  Z’0,  Z”v,  and  Kv  are  defined  by 
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Here  k  is  the  free  space  wave  number  and  kn,  is  the  transverse  wave  number  within  the 
waveguide.  The  equations  for  the  transverse  electric  and  magnetic  fields  beyond  the  aperture  (z  > 
0)  can  be  defined  via  a  Fourier  transform  as 
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where  the  ~  denotes  that  these  components  exist  in  the  spectral  domain.  In  order  to  determine 
E_T(kx,  ky)  and  Hr(ky,ky)  ,  we  start  by  equating  equations  5  and  7  at  z  =  0,  which  is  justified  by 

the  tangential  boundary  condition 

+00  +00 
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where  [i  =  x-x0  +  yy0.  Notice  that  on  the  right  side  of  equation  8  the  k  •  r  from  equation  7  has 
become  k*  g_aX  the  aperture  due  to  z  =  0.  The  following  orthogonality  equations  (5) 
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allow  us  to  simplify  equation  8  in  terms  of  the  reflection  coefficients  and  spectral  components  of 
the  E-  and  H-fields.  Here  S  is  the  dimensions  of  the  transverse  plane  of  the  rectangular 
waveguide.  By  limiting  the  bounds  of  integration  in  equation  9,  we  enforce  the  boundary 
condition  that  Ej(r)  =  0  on  the  surface  of  the  flange.  Using  the  following  general  substitution 
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we  can  rewrite  equation  8  as  the  following  system  of  equations: 
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This  set  of  equations  are  used  to  solve  for  ET(kx,ky)  and  Hr(kx,ky)  .  The  solutions  to  the 

surface  integrals  of  equation  10  can  be  determined  in  closed  form  and  are  calculated  in  the 
appendix.  In  order  to  get  equation  1 1  into  matrix  form,  we  need  to  add  like  ro  terms  together 
from  equation  1 1 .  We  can  now  rewrite  the  entire  system  of  equations  using  the  following  two 
expressions: 
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3.  Formulation  of  the  MDM  Method 


The  basis  of  the  MDM  method  is  the  system  of  equations  described  by  equation  12,  which 
populates  the  matrix.  By  representing  the  integrals  in  (kx,ky)  space  as  a  Riemann  sum  over  kx  and 
ky,  we  can  formulate  the  MDM  equation. 


3,1  Derivation  of  MDM  Equation 

To  simplify  equation  12,  we  can  write  equation  7b  in  terms  of  equation  7a  leaving  a  single 
unknown.  Starting  with  the  expression  (5), 
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we  can  rewrite  this  as  the  matrix  equation 
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Evaluating  the  cross  products  and  dot  product  of  equation  13  yields 


KAK’K)  = 


a)^o 


0 

k„  0 


KK 

Ko 

_K 

K 


Ko 

KK 

K„ 


E,(kx,kv) 


Ey(kx,kv ) 


1 


G>MoKo 


-KK  ~(K-K) 


k  -k 


KK 


Ex(kx,ky) 
Ey(kx,k  ) 


=  a.et(k,k) 


(14) 


Equation  12  represents  a  shorthand  notation  to  represent  the  system  of  equations  that  will 
construct  our  MDM  equation.  However,  to  properly  explain  how  these  equations  represent  a 
matrix  equation,  it  is  better  to  think  of  these  equations  in  longhand  notation: 
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This  set  of  equations  allows  us  to  solve  for  ET(kx,k  ) .  If  we  redefine  the  integrals  as  a  Riemann 
sum,  then  equation  15  becomes 
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(16b) 


M  '«*  ( K  >  K )  +  M  'oy  ( K » K )  Z0 


We  can  now  write  the  MDM  equation.  In  doing  so,  each  value  of  v  representing  a  different  mode 
inside  the  waveguide  represents  the  MDM  row  index.  Since  we  need  an  invertible  matrix  to 
obtain  a  solution,  v  also  represents  the  discrete  index  of  kx  and  kv  in  free  space,  which  becomes 
the  column  index  of  our  matrix.  This  creates  a  square  matrix  with  a  total  of  L  samples  of  kx  and 
kY  as  well  as  L  modes. 
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2  L 


2  L 


xo  -  direction 

v  -direction 
1—0 

M 

V(*n) 

M  \x(kn)  ••• 

M  1a  (kn ) 

M  1  y  (kn  ) 

M  ly  (kn )  •  •  • 

M  "l y  (krL  ) 

TE 

M' 

2a  (kn ) 

M  ix  (kn)  ■■■ 

M  "lx  ( kn  ) 

M  "ly  (kn  ) 

M  "ly  (kn )  •  •  • 

M  "ly  (krL  ) 

M' 

" Lx  (kn) 

M " Lx  (kn )  •  •  • 

M  " Lx  ( krL  ) 

M  Ly  (kn  ) 

M " Ly  (kn )  •  •  • 

M"Ly(kTL) 

M 

1a  (^ri  ) 

M  'lx(kn) 

M  1a  (kn. ) 

M  \y(kn) 

M  1  y  (kn)  ••• 

M  i y  (krc ) 

TM 

M 

2y{kn) 

M  2a  ( kn )  •  •  • 

M  ix  (kn ) 

M  ly  (kn) 

M  ly  (kn) 

M  'ly  (klL  ) 

M 

Nx  (  kn ) 

M  \x  (kn )  •  •  - 

M  ' Lx  (krL  ) 

M  Ly  (kn  ) 

M  Ly  (kn  )  •  •  • 

M  \y  (krL ) 

" 

Ex(k  n) 

2 

Ex  ) 

0 

Ex  ( krL ) 

• 

~ 

Eyikn) 

Eyikn ) 

- 

Ey(k  TL) 

0 

(17) 


In  this  notation,  v  has  taken  the  value  of  a  single  integer;  however,  each  row  index  truly 
represents  a  unique  wavenumber  (m,n)  pair.  The  MDM  of  equation  17  has  four  separate 
quadrants.  Quadrant  1  corresponds  to  the  matrix  elements  that  represent  the  TE  modes  inside  the 
waveguide  in  the  x-direction,  quadrant  2  corresponds  to  the  matrix  elements  that  represent  the 
TE  modes  inside  the  waveguide  in  the  y-direction,  quadrant  3  corresponds  to  the  matrix  elements 
that  represent  the  TM  modes  inside  the  waveguide  in  the  x-direction,  and  quadrant  4  corresponds 
to  the  matrix  elements  that  represent  the  TM  modes  inside  the  waveguide  in  the  y-direction 
where  (kj,,)  is  used  to  represent  the  (kx,ky)  pair.  The  reason  the  x-components  and  the  y- 
components  of  the  MDM  equation  have  to  be  separated  is  so  that  the  solutions  to  Ex(kx ,  ky  )  and 

Ey{kx,ky )  can  be  solved  for  individually.  Each  quadrant  is  L  x  L  in  dimension  yielding  a 

2 L  x  2 L  square  matrix.  The  solutions  to  Ex(k,kv)  and  Ey(kx , ky )  are  L  element  vectors.  Note 

that  the  right-hand  side  vector  of  equation  17  is  zero  except  for  the  first  element,  which 
corresponds  to  the  incident  TEio  mode  in  the  waveguide.  This  comes  directly  from  the  d„,v  in 
equation  12. 

As  with  any  discrete  representation  of  a  continuous  function,  L  must  be  large  enough  to  ensure 
an  accurate  representation  of  the  original  function.  However,  a  large  L  means  that  many  more 
modes  must  be  used  in  the  MDM  equation  than  are  truly  necessary  to  accurately  determine  the 
field  inside  the  waveguide.  This  can  lead  to  a  singular  matrix,  which  by  definition  is  not 
invertible.  Therefore,  in  solving  equation  17  we  use  singular  value  decomposition  (SVD)  to 
determine  the  pseudo-inverse  of  the  MDM  equation.  The  SVD  method  is  fully  described  in  most 
linear  algebra  texts  ( 6 ). 

3.2  Representation  of  kx  and  ky  When  z  >  0+ 

This  section  describes  how  to  represent  the  values  of  kx  and  ky  in  the  MDM  equation.  Since  it  is 
desirable  to  represent  the  stationary  phase  approximation  of  the  far  field  in  spherical  coordinates, 
we  must  map  kx  and  kv  to  (r,0,(j)).  The  stationary  phase  approximation  is  well  known  and  widely 
used  throughout  the  literature  (5,  7,  8).  The  equation  is  repeated  here  for  convenience: 
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(18) 


.  ke 


-jkr  r 


E(r,0,<l>)x  j  — -  9a  yEx{kx,ky^  cos  </>  +Ey  {kx,ky^  sin  (j) 


+  (f>  (^-Ex(jix,ky^sm(j)  + Ey(kx,ky^cos(j)^cos6 


Mapping  kx  and  ky  to  spherical  coordinates  yields 

kx  =  ka  sin(#)cos(0) 

ky  =  k()  sin (6>) sin (0) 


(19a) 

(19b) 


In  order  to  get  a  hemisphere  mapping  of  the  radiated  far  field  in  the  propagation  direction,  we  are 
interested  in  values  of  -n/2  <  9  <  n/2  and  0  <  ([>  <  n/2.  After  substituting  these  values  of  9  and  <f> 
into  equation  19,  we  get  a  trajectory  of  kx  and  kv  onto  a  circle  of  radius  k0,  as  shown  in  figure  2. 


Figure  2.  Plot  of  the  values  of  kx  and  kr  obtained  for  <j)  =  0. 

Figure  2  shows  all  the  kx  and  ky  values  obtained  for  </)  =  0,  -n/2  <  6  <  n/2,  and  A 0=  ML,  where  L 
corresponds  to  the  size  of  each  quadrant  in  the  MDM  equation.  The  angle  of  (j)  is  represented  in 
figure  2  as  the  angle  between  the  kx  and  ky  axes.  Since  <j>  =  0,  all  the  kx  and  kv  values  fall  on  the 
ky  =  0  axis.  If  we  use  <j>  =  rc/4  to  calculate  kx  and  ky  as  in  figure  3,  then  we  see  that  the  values  of  kx 
and  kY  fall  on  a  trajectory  that  makes  an  angle  of  jt/4  with  the  ky  =  0  axis. 

Using  this  technique,  we  can  generate  the  radiated  far  field  for  any  value  of  0  <  <f>  <  n/2.  Note  that 
in  both  figures  2  and  3  an  equal  spacing  between  values  of  9  does  not  result  in  an  equal  spacing 
in  kx  and  ky.  Also,  any  value  of  kx  and  ky  that  falls  on  the  k()  radius  yields  a  value  of  \kj\  =  k0, 
which  corresponds  to  K  =  0.  Any  value  of  kx  and  kv  that  falls  beyond  the  k0  radius  corresponds  to 
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an  imaginary  value  of  K.  These  values  represent  attenuating  modes  and  are  not  used  in  the 
calculations  for  the  incident  TEio  mode. 


Figure  3.  Plot  of  the  values  of  kx  and  kr  obtained  for  <j>  =  jt/4. 


4.  Calculation  of  Far  Field  and  Comparison  to  Simulation 


This  section  compares  the  far- field  radiation  pattern  calculated  from  the  MDM  method  described 
in  section  3  to  those  generated  using  CST  Studio  Suite  2012.  The  MDM  calculations  were 
performed  using  Matlab. 

4.1  Description  of  CST  Model 

Figure  4  depicts  the  model  used  to  simulate  the  semi-infinite  rectangular  waveguide  with  an 
infinite  flange.  Figure  4a  gives  the  transverse  dimensions  of  the  waveguide,  where  a  =  a/2  at 
200  MFIz  and  b  =  a/2.25.  As  long  as  a  >  b,  the  dominant  mode  will  be  the  TEio  mode.  The  cutoff 
frequency  for  the  propagation  of  the  dominant  mode  is  200  MHz  and  the  cutoff  frequency  for  the 
next  mode  to  propagate  is  400  MHz.  For  frequencies  below  200  MHz,  no  modes  will  propagate 
in  the  waveguide,  and  for  frequencies  above  400  MHz,  more  than  one  mode  will  propagate  in  the 
waveguide.  Since  the  method  described  in  section  2  corresponds  to  a  waveguide  that  is  infinite  in 
one  direction  and  ends  at  a  radiating  flange  at  z  =  0,  placing  a  waveguide  port  at  the  end  of  the 
waveguide  mimics  this  setup,  as  shown  in  figure  4b.  The  waveguide  port  will  absorb  any 
reflections  of  additional  modes  from  the  aperture  to  ensure  that  only  the  propagating  TEio  mode 
will  contribute  to  the  far  field  radiation.  The  direction  of  the  arrow  in  figure  4b  shows  the 
direction  of  propagation  for  the  TEio  mode. 
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Figure  4.  (a)  Transverse  plane  of  the  waveguide  aperture.  The  waveguide  has  dimensions  |a|  x  |b| 
and  is  surrounded  by  an  infinite  conducting  flange,  (b)  This  shows  the  z-direction  of  the 
waveguide  with  the  waveguide  port  included. 


Figure  5  shows  the  distribution  of  the  mode  generated  by  the  waveguide  port  as  a  cosine 
distribution  across  the  long  dimension  of  the  waveguide  aperture  and  has  units  of  volts/meter 
(V/m).  This  is  the  expected  mode  distribution  for  the  TEio  mode  (P).  Notice  that  the  mode 
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distribution  peaks  and  is  symmetric  about  x  =  0  and  y  =  0  as  expected.  The  mode  distribution 
does  not  vary  in  the  z-direction  because  the  mode  is  propagating  and  not  attenuating. 


Figure  5.  Illustration  of  the  mode  distribution  in  the  rectangular 
waveguide  of  the  infinite  flange  simulation  in  V/m. 


4.2  Result  Analysis 

The  calculations  of  the  far- field  radiation  patterns  are  determined  for  a  frequency  of  300  MHz 
This  frequency  was  chosen  because  it  stands  farthest  away  from  both  the  cutoff  that  will  not 
allow  the  TEio  mode  to  propagate  and  the  cutoff  that  will  allow  multiple  modes  to  propagate. 
The  far- field  radiation  patterns  in  the  90  and  (ji,,  directions  are  calculated  from  equation  18  as 


Ex  (kx ,  ky ,  z  =  O)  cos  (f)  +  Ey  [kx  ,ky,z  =  0)  sin  </) 

+  t-o  (“•/*  sin  0  +  fy  COS  (f)  COS  6 


(20) 


Figure  6a  shows  the  patterns  of  the  far  field  in  the  Fq  and  F(j,  directions  plotted  on  a  polar  graph, 
whereas  figure  6b  shows  the  patterns  of  the  far  field  in  the  Fg  and  directions  plotted  linearly. 
These  plots  assume  sr  =  1  (relative  permittivity)  and  /ir  =  1  (relative  permeability)  inside  the 
waveguide.  We  can  see  from  inspection  that  the  results  of  solving  the  MDM  equation  agree  very 
closely  with  the  results  generated  by  the  simulation  using  CST  Studio  Suite. 
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Figure  6.  (a)  Polar  plot  (b)  Linear  plot  -  of  the  far-field  Fg  and  /y 
normalized  radiation  patterns. 


One  thing  to  note  is  that  when  using  the  SVD  method  in  these  computations  the  number  of 
singular  values  used  to  generate  the  pseudo-inverse  of  the  MDM  plays  a  crucial  role.  A  matrix 
with  dimensions  2 L  x  2 L  will  have  2 L  singular  values.  Many  of  the  singular  values  will  have 
magnitudes  approaching  zero.  These  values  should  not  be  used  or  they  will  affect  the  accuracy  of 
the  numerical  results.  On  the  other  hand,  if  one  has  multiple  singular  values  with  useable 
magnitudes,  then  eliminating  any  of  them  from  the  calculations  will  also  affect  the  results. 

Figure  7  shows  a  plot  of  the  singular  values  in  descending  order  for  the  MDM  calculation 
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resulting  in  the  patterns  of  figure  6.  The  number  of  singular  values  used  for  this  calculation  is  4. 
Generally  only  the  first  few  singular  values  are  needed. 


Index 


Figure  7.  Plot  of  the  singular  values  of  the  MDM  in  descending  order. 


5.  Conclusions 


Existing  methods  for  analyzing  the  radiating  infinite  flange  require  the  computation  of  the  fields 
at  the  aperture.  These  fields  then  have  to  be  transformed  to  the  spectral  domain  before  any  far- 
field  calculations  can  be  made.  This  report  derives  a  new  approach  called  the  MDM  method  that 
allows  for  the  direct  computation  of  the  radiated  far  fields.  The  result  is  a  matrix  equation  that 
directly  solves  for  the  spectral  components  needed  for  the  far-field  stationary  phase 
approximation.  The  results  of  the  MDM  method  were  successfully  compared  to  the  far-field 
radiation  patterns  of  an  infinite  flange  generated  by  commercial  numerical  modeling  software  for 
the  same  problem. 
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Appendix.  Derivation  of  Surface  Integrals  for  MDM  Equation 


This  section  focuses  on  the  solutions  to  the  surface  integrals  in  equations  10-12.  If  we  assume 
that  the  dimensions  of  our  radiating  aperture  are  known  and  that  the  dimensions  are  rectangular 
with  values  of  a  >  b,  respectively,  then  these  surface  integrals  can  be  solved  in  closed  form.  The 
solutions  yield  constants,  which  are  used  to  populate  the  MDM  in  equation  17. 

A-l  Surface  Integrals  for  TM  Modes 

We  first  solve  the  surface  integral  that  arises  from  the  electric  field  generated  by  the  TM  modes 
in  a  rectangular  waveguide.  The  solutions  for  the  incident  and  non-incident  modes  have  the  same 
form,  so  we  use  the  symbol  v  to  denote  either  N  or  M. 

We  start  with  the  following  equation,  which  represents  the  two  dimensional  surface  integral 
corresponding  to  the  TM  mode  vector  for  the  E-field, 

b  a 

e'  (kx,ky)=  |  I  e  ,tr'/’ e'  (x,y)dxdy  (A-l) 

_b  a 
2  2 


where_e  ’„(x,y)  is  defined  by  equation  2a.  We  can  extract  the  x-component  and  solve  for  it 
individually: 


(A-3a) 


a 


2 


(A-3b) 


then  equation  A-3a  has  a  known  closed-form  solution  when  ky  ±  0  and  n  ±  0  from  (10): 


17 


w=- 


-jkyy 


^  qn^ 


~k>'+'  b 


-jkv  sin 


^  nn  (  b\\  b 

—  v  +  — 

V  b  V  2  JJ 


nn 

\bj 


cos 


b nn  f 


b 

y  +  - 

V  u  V  2  j) 


V  o  j 


k  nn^ 


\  o  j 


K  - 


r  nn  A 

\b  J 


b 

-jK- 


(-1  r- 


+jK 


k  ^  0,  n  =  1, 2, 3 . .. 


When  ky  =  0  and  n±  0,  then 

b 

/,<*,)=  J 


b nn  b 


sin 


b 

y  +  - 

V  y  V  2yy 


dy 


i-(- 

177T  V 


kv  =0,  n  =  1,2,... 


Finally,  when  n  =  0  then  Ii(ky)  =  0.  Similarly,  when  kx  #  0  and  m  ^  0,  then 
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and  when  /«  =  0  and  kx  ±  0,  then 
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Ii(kx)=  J  e~jk*xdx  =  -jk 
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Finally,  when  kx  =  0  and  m  =  0,  then 
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~C3  |  <N 


z 

I2(kx )  =  J  dx  =  a  kx  =  0,  m  =  0 


Therefore,  the  final  form  of  equation  A-3  is 
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A:  ^  0,  m  =  1,2,3... 


A\.  =  0,  m  =  1,2,3... 


kx  *  0,  /«  =  0 


A  =  0,  /«  =  0 


(A-4a) 


(A-4b) 


Now  that  we  have  closed- form  solutions  for  the  integrals  /,  and  /2 ,  we  can  determine  the  final 
form  of  A-2  as 


e'(kx,kv)  =  -A'u—Il(k)I2(kx) 


(A- 5) 


Similarly,  we  can  extract  the  y-component  from  equation  A-l  and  solve  for  it  individually  as 


■  ..  ,  .  rm  r  -ikX  .  ( m7r  (  ,  f  -jkv  (nn(  b\]  ,  , . 

?  (kr,k„)  =  -A  n —  e  sin  -  x  +  —  dx  e  y'  cos  —  v  +  —  dv  (A-6) 

v  h  '  n  ?  J  /,  1  -  9  v  7 

D  a  v  a  V  b  v  0  V 
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(A-7a) 
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(A- 7b) 
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then  using  the  same  analysis  we  used  to  solve  the  integrals  of  equation  A-3  we  can  write 
equation  A-7  as 
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777  =  0 

kx  =  0, 777  =  1,2,3,.. 


(A- 8  a) 


(A-8b) 


Now  that  we  have  closed-form  solutions  for  the  integrals  I3  and  I4,  we  can  determine  the  final 
form  of  A-6  as 


e'  (kx,ky)  =  -A'v^I3(ky)I4(kx) 
y  b  (A-9) 

Now  we  solve  the  surface  integral  that  arises  from  the  magnetic  field  generated  by  the  TM 
modes  in  a  rectangular  waveguide.  We  begin  with  the  following  equation,  which  represents  the 
two-dimensional  surface  integral  corresponding  to  the  TM  mode  vector  for  the  H-field, 

b  a 
2  2 

li 1  ( kx,ky)=  |  Je  Jb’Pfoi  (x,y)dxdy  (A- 10) 

b  a 
2  2 

where  h  ’„(x,y)  is  defined  by  equation  2b.  We  can  extract  the  x-component  and  solve  for  it 
individually: 


20 


h'Akx’ky)  =  A 


H7T 


2 

f 

i 

i 

** 

* 

S' 

7  niK 

( 

a* 

_ i 

2 

J 

e  Jk>’y  cos 

7  n7T 

<3  1  < 

+ 

\ 

dy 

J 

a 

1 

l  a 

l  2 ) 

j 

J 

b 

_ 

l  b 

f  2/ 

) 

2  2  (A-ll) 

We  see  that  these  are  the  same  integrals  as  those  of  equation  A-6,  and  we  can  write  the  solution 
to  A- 1 1  as 

h'Jkx,ky)  =  A^I3(ky)I4(kx) 

b  (A- 12) 

Shnilarly,  we  can  extract  the  y-component  from  equation  2b  and  solve  for  it  individually: 
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(A- 13) 


We  see  that  these  are  the  same  integrals  as  those  of  equation  A-2,  and  we  can  write  the  solution 
to  A-13  as 


nm 


h „  (K » K  )  =  ~A'v - !\(kv  )72  (kx  ) 


(A- 14) 


We  can  now  rewrite  equation  12b  in  terms  of  equations  A-5,  A-9,  A- 12,  and  A- 14  as 
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(A- 15) 


Equation  A- 15  represents  the  final  form  that  we  need  to  determine  our  system  of  equations  for 
the  TM  modes. 


A-2  Surface  Integrals  for  TE  Modes 

Now  we  solve  the  surface  integral  arising  from  the  electric  field  generated  by  the  TE  modes  in 
the  rectangular  waveguide.  We  begin  with  the  following  equation,  which  represents  the  two- 
dimensional  surface  integral  corresponding  to  the  TE  mode  vector  for  the  E-field, 

b  a 

e "  (kx ,ky)  =  J  }  J~T e (x, y)dxdy  (A- 1 6) 

b  a 
2  2 

where  e  ’  ’v(kx,ky) is  defined  by  equation  2c.  We  can  extract  the  x-component  and  solve  for  it 
individually: 
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(A- 17) 


We  see  that  these  are  the  same  integrals  as  those  of  equation  A-2  and  we  can  write  the  solution  to 
A-17  as 


« V  (K  X,,)= 2j=b  (krVi  (K ) 

bylab 

Similarly,  we  can  extract  the  y-component  from  equation  2c  and  solve  for  it  individually: 


(A- 18) 
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(A- 19) 


We  see  that  these  are  the  same  integrals  as  those  of  equation  A-6  and  we  can  write  the  solution  to 
A-19  as 


^  /  11  1  7T 

e\(kx,ky)  =  --^I3(ky}IA(kx) 

a^ab  (A-20) 

Now  we  can  solve  the  surface  integral  arising  from  the  magnetic  field  generated  by  the  TE 
modes  in  the  rectangular  waveguide.  Beginning  with  the  following  equation,  which  represents 
the  two  dimensional  surface  integral  corresponding  to  the  TE  mode  vector  for  the  H-field, 

b  a 
2  2 

h”  ( kx,ky)=  J  Je  (x,y)dxdy  (A-21) 

b  a 
2  2 


where  _h”v(kx,kv)  is  defined  by  equation  2d.  We  can  extract  the  x-component  and  solve  for  it 
individually: 
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We  see  that  these  are  the  same  integrals  as  those  of  equation  A-6,  and  we  can  write  the  solution 
to  A-22  as 


nui 


h\SK,ky)  =  A\—I,(kv)IA(kx) 


3V'-yU4VvxJ 


(A-2  3) 
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Similarly,  we  can  extract  the  y-component  from  equation  2d  and  solve  for  it  individually: 

a  b 


2  2  (A-24) 

We  see  that  these  are  the  same  integrals  as  those  of  equation  A-2,  and  we  can  write  the  solution 
to  A-24  as 


h  \y  ’ky)  =  A  —  /j  (ky  )I2  (kx  ) 


(A-25) 


We  can  now  rewrite  equation  12a  in  terms  of  equations  A- 18,  A-20,  A-23,  and  A-25  as 


e"Jkx,ky)Xo+e"  Jkx,ky)yo  ,  ET(kT)J 

+Z  "u  [Aj1  ’  ky  )E  +  h  "ljy  (K  ’  ky  )y0  ]  ’  ET(kT) 


(A-2  6) 


Equation  A-26  represents  the  final  form  that  we  need  to  determine  our  system  of  equations  for 
the  TE  modes. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


FDTD 

finite-difference  time-domain 

MDM 

modal  decomposition  matrix 

MHz 

megahertz 

SVD 

singular  value  decomposition 

TE10 

transverse  electric  with  mode  number  (1,0) 

TM 

transverse  magnetic 

V/m 

volts/meter 
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